Mon. Not. R. Astron. Soc. 000,[T][T0](2007) Printed 4 February 2008 (MN WpX. style file v2.2) 



Selection of ThAr lines for wavelength calibration of echelle spectra 
and implications for variations in the fine-structure constant 
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ABSTRACT 

Echelle spectrographs currently provide some of the most precise and detailed spectra in as- 
tronomy, the interpretation of which sometimes depends on the wavelength calibration accu- 
racy. In some applications, such as constraining cosmological variations in the fundamental 
constants from quasar absorption lines, the wavelength calibration is crucial. Here we detail an 
algorithm for selecting thorium-argon (ThAr) emission lines for wavelength calibration which 
incorporates the properties of both a new laboratory wavelength list and the spectrograph of 
interest. We apply the algorithm to the Very Large Telescope Ultraviolet and Visual Echelle 
Spectrograph (UVES) and demonstrate a factor of >3 improvement in the wavelength cali- 
bration residuals (i.e. random errors) alone. It is also found that UVES spectra calibrated us- 
ing a previous, widely distributed line-list contain systematic ±30-75 ms~' distortions of the 
wavelength scale over both short and long wavelength ranges. These distortions have impor- 
tant implications for current UVES constraints on cosmological variations in the fine-structure 
constant. The induced systematic errors are most severe for Mg/Fe n quasar absorbers in the 
redshift range 1.2 <z a bs S>2.3, with individual absorbers studied by recent authors containing 
systematic errors up to 4 times larger than quoted statistical errors. 
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1 INTRODUCTION 

Echelle spectrographs are the most common instrument for record- 
ing high-dispersion astronomical spectra. Since many echelle 
diffraction orders can be cross-dispersed and thereby recorded si- 
multaneously on rectangular-format media [e.g. modern charge- 
coupled devices (CCDs)], more-or-less continuous wavelength 
coverage over much of the optical range can be achieved whilst 
maintaining high resolving powefl To fully exploit the veloc- 
ity precision available in high-resolution spectra, precise wave- 
length calibration over much of the optical range is required. 
This is now usually done by comparison with exposures of a 
thorium-argon (ThAr) hollow-cat hode emission-line lamp (e.g. 
B rec kinridge. Pierce & Stolll 19751) . 

After iMeggerj dl957T) first suggested that Th transitions 
might act as reliable wavelength standards, many authors (e.g. 
lGiacchettietai]|l970h measured Th laboratory wavelengths and 
confirmed their suitability for such purposes. There are several rea- 
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1 Resolving power R = A /FWHM, where FWHM is the instrumental pro- 
file's full-width at half-maximum. 



sons the ThAr spectrum has proven an enduring standard for as- 
tronomical spectroscopy: (i) Most optical ThAr transitions have 
measured laboratory wavelengths, most of which are resistant to 
pressure and Stark shifts; (ii) Th has one stabl e isotope and ter- 
restri al Ar is almost entirely composed of 40 Ar dRosman & Taylor] 
1998). Also, both Th and 40 Ar have even nuclei. Therefore, no un- 
resolved isotopic or hyperfine structure in the ThAr lines is ex- 
pected; (iii) Most ThAr lines are unresolved at most practical re- 
solving powers, i.e. R < 150000; (iv) The ThAr line-density is 
high, > 1 line per A, over most of the optical range; (v) ThAr lamps 
are fairly cheap, are commercially available and have lifetimes of 
typically > 1 yr. Other calibration sources include iodine absorption 
cells (e.g. iMarcv & Butledfl992r) but these degrade the object sig- 
nal and can only be used over a comparatively short wavelength in- 
terval (typically ~5000-6500 A). Another, potentially ultra-stable, 
absolute calibration system might be based on t he new technology 
of laser frequency combs dMurphv et alj|20073) . a possibility now 
under development. 

Despite the distinct advantages, there are, however, sev- 
eral drawbacks to using ThAr spectra for wavelength calibration. 
Firstly, different ThAr lines have widely different intensities, far 
beyond the typical dynamic range of modem CCDs (i.e. a factor of 
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~50000). Secondly, the wavelength distribution of lines is highly 
non-uniform; it turns out that most ThAr lines are blended with 
relatively strong nearby lines, even at R < 150000. It is therefore 
important to carefully select which ThAr lines are to be used when 
calibrating echelle spectra. This task is complicated by the fact that 
a given line's intensity (relative to other lines) varies from lamp 
to lamp and even within the one lamp over time. ThAr lines arise 
from a variety of Th and Ar ionization states and so, as the lamp 
ages and/or the operating current and pressure are varied, the rela- 
tive strengths of different lines can change. Thus, ideally, detailed 
selection of ThAr lines for wavelength calibration is best carried 
out separately for each ThAr exposure. However, in practice, most 
spectrograph-specific data reduction pipelines carry out only sim- 
ple and minimal line selection, thereby assuming that a detailed 
pre-selection of lines has already been carried out. That is, for most 
practical purposes, it suffices to have a ThAr line-list pre-selected 
for the spectrograph of interest. 

In this paper we detail an algorithm to pre-select ThAr lines 
for any echelle spectrograph and we apply the algorithm to the Ul- 
traviolet and Visual Echelle Spectrograph (UVES) on the European 
Southern Observatory's (ESO's) Very Large Telescope (VLT) on 
Cerro Paranal, Chile. ESO have already pre-selected ThAr lines to 
calibrate UVES and the line-list is included with the ESO pipeline 
fo r UVES data reduction. T h eir se lection der ives from the results 
of iDe Cuyper & Hensbergi d 19981 hereafter DCH98) who con- 
structed simu lated ThA r emission spectra based on line tables in 
the literature. IX'l 198 report 'effective' wavelengths for selected 
ThAr lines given the expected degree of line blending based on 
the relative line intensities in the literature. The effective wave- 
lengths therefore depend on resolution; the ESO list is based on 
the R = 1000 00 results of|DCH983 However, the ESO list is not 
identical to the DCH98 list and seems to include > 500 lines which 
were rejected by DCH98 because no effective wavelengths could 
be found. Furthermore, all wavelengths are only given to 3 decimal 
places (in units of A) and, in many cases, the wavelengths were 
truncated, rather than rounded, from 4 decimal places. 

Such treatment of the ThAr line-list can cause non-negligible 
systematic effects in several applications. In this paper we fo- 
cus on the effects it has had on current constraints on cosmo- 
logical variations in the fine-structure constant, ar em , derived from 
VLT/UVES quasar absorption spectra. If a cm were different in 
the high-redshift absorption clouds, the relative wavelengths of 
different metal-line transitions would differ from those observed 
in laboratories on Earth. This effect provides a precise observa- 
tional probe of possible variations in or e m over cosmological time- 
and distance-scales called the 'ma ny-multiplet' (MM) method 
jDzuba et alJl 19991 ; I Webb et alJI 19991) . Analyses of a large sample 
of quasar spectra from the Keck High Resolution Echelle Spec- 
trograph (HI RES) have so far provided the first evidence for vari- 
ation in a em jWebb et alj|l999l 1200 ll ; iMurphv et all 1200 ll , l2003h 
with the most recent results indicating a smaller a m in the absorp- 
tion clouds at the fra ctional level of ha /a = (-0.57 ± 0.11) x 10~ 5 
(Mur phy etal]|2004l) . Clearly, such a potentially fundamental re- 
sult must be refuted or confirmed with many different telescopes 
and spectrograph s. Some first attempts with VLT/UVES have now 
been made (e.g. IChand e"tai]|2004 l20o3 ; iLevshakov et alj|200^) . 
generally yielding no variation in a em . While detailed analyses of 
HIRE S ThAr spectra have show n no signs of strong systematic er- 
rors ( Mur phy etal.ll200ll , l2003l) . no similar general assessment has 
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yet been made for the bulk of the UVES constraints on ha /a. Given 
the errors in the ESO ThAr line-list, it is imperative to do so. 

This paper is organised as follows. In Section [2] the various 
ThAr transition tables in the literature are combined according to 
their different strengths and weaknesses. Section [3] details the line 
selection algorithm. Section|4]describes the final ThAr line-list for 
use in calibrating VLT/UVES spectra and shows the distortions of 
the wavelength scale introduced by the errors in the ESO line-list. 
Section [5] assesses the general implications these distortions have 
for current UVES constraints on a varying fine- structure constant. 



2 SYNTHESIS OF EXISTING TH AND AR LINE-LISTS 

Before selecting which ThAr lines are to be used to calibrate 
UVES, we first require a list of laboratory wavelengths for as many 
features as possible which appear in measured UVES ThAr lamp 
spectra. Most features, especially below 7000 A, are either known 
to be due to Th or are unidentified (but probably due to Th), while 
~7 per cent of the features are from Ar and < 1 per cent are from 
'contaminant' species such as Nai, Mgi, Can and Fei. The Th and 
Ar lines, and even some of the contaminant lines, are catalogued 
in various atlases derived from painstaking (but usually rather old) 
laboratory work. These atlases provide our starting point. However, 
every ThAr lamp gives a somewhat different spectrum and there al- 
ways appear additional lines which cannot be found in any atlas, 
even as 'unidentified' lines. These are probably due to additional 
contaminant ions and molecular species in the lamp. This means 
that our prior knowledge of the ThAr spectrum from any given 
lamp is, at best, incomplete, thereby necessitating the selection pro- 
cedure detailed in the next section. 

There is no single ThAr atlas covering the whole wavelength 
range of UVES, -3030-10540 A. How ever, there is one atla s of Th 
lines onl y whic h doe s cover this range, Palm er & Englemanl Jl983l. 
hereafter |PE83J). The |PE83j absolute velocity preci sion va ries from 
~15 to ~120ms~' depending on the line intensity. |PE83l identify a 
large number of Th i, n and m lines in their spectrum but a large 
number of unidentified lines remain, most of which are probably 
due to Th. Furthermore, on comparison with real ThAr lamp spec- 
tra, one notices that there are still several thousand lines which can- 
not be accounted fo r by Ar or cont a minan t species. 

Very recently, lLovis & Pepej d2007l. hereafter ILP07I) used a 
large library of ThAr spectra taken over one month with the High 
Accuracy Radial Velocity Planet Searcher (HARPS) spectrograph 
on the 3.6-m ESO telescope to improve the measurement precision 
of many ThAr lines, at least over the HARPS wavelength ra nge of 
3785-6915 A. They identified in their spectra lines from the |PE83l 
atlas which showed no strong wavelength variations with time and 
obtained line positions with a root-mean-square variation (RMS) of 
typically 5-10 m s -1 . Thus, while bootstrapping their overall wave- 
length scale to that of PE83, they were able to correct the wave- 
lengths of individual PE83 lines, especially weaker lines. More- 
over, they were able to measure the wavelengths of all other fea- 
tures in their ThAr spectra, again at the 5-10 ms~' precision level. 
They removed from their list lines which saturated their CCD, (very 
few) lines which appeared in the |PE83l catalogue but which were 
not detected in the HARPS spectra and, most importantly, they re- 
moved lines which were either closely blended with other lines at 
the HARPS resolving power of R = 120000 or were observed to 
change position with time. The latter indicates either that the line 
experiences significant shifts with changing lamp pressure or cur- 
rent or that the line is actually a blend and that the relative inten- 
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sities of the blended lines varies with changing lamp conditions. 
In this work we use all lines from the lLP07l atlas excluding those 
which have relative wavelength errors amou nting to > 120 m s ; a 
total of 8 219 ou t of 8442 lines are used from lLPOl 

The lLP07l catalogue ranges from 3785 t o 6915 A. In selecting 
lines for use in calibrating UVES, we use the PE83 line list for Th 
and contaminant lines outside this range. However, it is important 
to begin with a list which includes the known Ar lines as well so 
that (i) the maximum number of known blends can be identified and 
removed and (ii) the maximum number of lines are available for 
calibrating UVES. Three A r line-lists exist w ith similar wav e length 
precision - Ikorlenl jl973l hereafter |N73|) . IWhaling et all Jl995l 
hereafter Iw95h . and IWhaling et al.1 J2002l . hereafter fWpll) - ea ch 
of which has different advantages and disadva ntages . The lN73l at- 
las contains both Ari and Arn lines whereas |W95| contains Arn 
lin es on ly and lW02| contai ns An lin es only. Fewe r lines are listed 
by |N73l compared to lW95l an d lW02l since the lN73l experimen t was 
less sensitive to weak lines. The velocity precision achieved by N73 
for th e Ar i lines is better than that of W02 but the velocity precision 
ofthe lW95l Ar ii lines is better than N73's. Finally, there is a calibra- 
tion differen ce such that the Wha ling et alj wavenu mbers are large r 
than lNorlerf s by a factor of S NW = 1 + 6.8 x 10~ 8 . IWhaling et all 's 
calibration should be more reliable since it uses calibration lines in 
the same rec orded spectra as the Ar n lines. Therefore, we use the 
IWhaling et al J calibration scale in synthesising the Ar line-lists. 

The different lists of An and Am lines are combined in the 
following ways. We do not consider the Arm lines from lW95l since 
there are very few of them and since their wavelengths have larger 
uncertainties (W. Whaling, private communication). 

• An: Following the recommendation in |W02| we used |N73| 
wavenumbers (scal ed by Snw) when available, othe rwiselW02l val- 
ues were used. For N73 lines, the intensity scale of lN73l was used 
as re- cast by |DCH 98 onto a logarithmic scale. For lW02l lines, the 
|W02| intensity scale was used. 

• Ar ii: |W95I wavelengths were used here since |W95| claim that 
they are relatively ins ensiti ve to pressure-shifts and that th eir values 
are more prec ise tha n lN73l 's. In t he fe w cases where lN73l reports an 
Arn line that W95 doesn't, the wavelength (scaled by Snw) 
was used. For lW95| lin es the lW95l intensit y scale was used (this is 
the same as the | W02| intensity scale). For |N73| values, the inten- 
sity scale of lN73| was used as re-cast by DCH98 onto a loga r ithmic 
scale. Note that lN73l uses the intensity scale of lMinnhager] d 19731) 
for Am lines above air wavelengths of = 7600 A and below 
'iair = 3400 A. This subtlety is neglected in the considerations be- 
low. 

After combining the iLPOl |PE83l lN73l |W95| and lwoj lists in 
the above way we obtain the list hereafter referred to as LPWN. 
Note from the above that the intensity scales for each species (e.g. 
Th i, Ar i, Ar n) are different and that the relative intensities of lines 
from different species depend on many factors such as lamp pres- 
sure and current and the relative partial pressures of Th and Ar gas 
used. 



3 THAR LINE SELECTION 

One should not simply use the above line-list, LPWN, to cali- 
brate UVES spectra because of the potent i al for line-blending. The 
Fourier transform spectra of |PE83l lN73l |W95| and |W02| and the 
HARPS spectra of lLP07l all have resolving powers R > 120000 
whereas the UVES resolving power is typically 35000-70000 and 



possibly as high as 110000. Therefore, a major part of the line se- 
lection that follows is the rejection of 'close' blends. The concept 
of 'close' clearly depends on the relative intensities of the blended 
lines and so it is necessary to put all lines in LPWN list onto the 
one intensity scale typical of UVES ThAr spectra. Furthermore, 
since there are always additional unknown lines in measured ThAr 
spectra, a given ThAr line can only be used for calibration if it is 
measured to be reliable in the UVES ThAr spectrum. For these rea- 
sons we have constructed a UVES ThAr spectrum for use in the line 
selection process. 



3.1 UVES ThAr spectrum 

ThAr spectra taken as part of regular UVES maintenance ac- 
tivities with a 076 slit and no CCD rebinning were retrieved 
from the ESO Science Archive. Exposures from several differ- 
ent standard UVES wavelength settings^ (346, 437, 580, 600, 760 
and 860-nm) were used to provide complete wavelength cover- 
age, with the exception of three small echelle order gaps red- 
wards of 10000 A (10084.43-10085.93 A, 10252.42-10256.95 A 
and 10426.10-10433.76 A). The UVES pipelintQ was used to ex- 
tract and wavelength calibrate the data. Several modifications to 
the pipeline, which will be reported elsewhere, were made to im- 
prove object extraction and wavelength calibration. The data were 
re-dispersed to a log-linear wavelength scale using uves POPLEF0, 
code specifically written to re-disperse and combine UVES data 
from multiple wavelength settings. Overlapping regions of spec- 
tra were cut away so that only one raw exposure contributed to 
the final spectrum at any given wavelength. No effort was made 
to flux-calibrate the final spectrum and the blaze-function of the 
echelle grating was still evident in the combined spectrum. How- 
ever, our results from the intensity re-scaling in the next section 
demonstrate that this is a minor consideration. The final spectrum 
has a resolving power of R as 70000 and the log-linear dispersion 
is set to 1.75 km s~ . 



3.2 ThAr line selection algorithm 

The above line-list, LPWN, is treated as the input catalogue to 
the following algorithm for selecting a final list for calibration of 
UVES. We make three passes through the algorithm, refining at 
each pass the line list used to calibrate the UVES ThAr spectrum 
and other parameters specified below. Figures pertain to the 
final pass through the algorithm. 

(i) Wavelength calibrating the UVES ThAr spectrum: The pro- 
cedure in Section [3~T1 for constructing the UVES ThAr spectrum 
requires an input line-list for the wavelength calibration stage. At 
the first pass we use the line-list provided by ESO with the UVES 
pipeline. At the second and third passes the output from the previ- 
ous pass through the algorithm is used. Note that the line-list used 
at the first pass is irrelevant; only an approximate wavelength scale 
needs to be established initially and this is greatly refined by sub- 
sequent passes through the algorithm. 

(ii) Gaussian fitting: Each line in the input list LPWN is 
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sures were taken on the 1 1th February 2004. 

4 http://www.eso.org/projects/esomidas 

5 http://www.ast.cam.ac.uk/~mim/UVES_popler.html 
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Figure 1. Examples of the inten sity sc aling procedure. Measured and catalogued intensities, / mca s and hin respectively, for unblended lines from each category 
- Th i from lLP07l and Th n from |PE83l in this figure - are compared to derive the scale-factor, a. All lines from that category are then scaled by a to the UVES 
intensity scale. 



searched for in the UVES ThAr spectrum and fitted with a Gaus- 
sian. The fit includes the 13 pixels centred on the pixel with maxi- 
mum intensity. A first-guess continuum level is defined by averag- 
ing the first and last 2 pixels in the window. A first-guess contin- 
uum slope is defined from the difference between these averages. 
A first-guess amplitude is defined as the maximum intensity minus 
the continuum level. A first-guess width is defined from the veloc- 
ity difference between the first and last pixels which have intensi- 
ties less than the continuum plus half the first-guess amplitude. The 
first-guess central velocity is calculated as the intensity-weighted 
velocity of the 3 central pixels. A 5-parameter Gaussian fit is per- 
formed and the best-fit parameters are used in the subsequent steps 
of the line selection algorithm. The initial guesses and Gaussian fit- 
ting procedure are identical to those we used in our modifications 
to the UVES pipeline. 

(iii) Intensity re-scaling: The different Th and Ar atlases com- 
bined above have different intensity scales. Furthermore, one would 
not expect to find the same relative intensities for lines of different 
ionization stages (e.g. An and Am) in the laboratory and UVES 
spectra since the ThAr lamps will have had different operating con- 
ditions (e.g. pressure, current etc.). We therefore aim to place all 
lines from the different ionic species (i.e. Thi/n/m, Ari/n, contam- 
inan t and unidentified spec ies) from e ach different intensity scale 
(i.e. iLPOl IPE83L lN73l and |W95|[W02|) on a single intensity scale 
directly related to the UVES ThAr spectrum. First, any pairs of 
lines within 13 km/s of each other are removed from the input line 
list LP WN. For each 'category' of line - for exa mple, Thi from 
ILPOl Thi from |PE83l Arn from lLP07l Am from lW95l etc. - the 
catalogued intensities are compared with those measured from the 
Gaussian fits performed on the UVES ThAr spectrum. The me- 
dian ratio of measured and catalogued intensity, a, is defined as 
the scale-factor. All lines from this category are then scaled to the 
UVES intensity scale by multiplying their catalogued intensities by 
a. Two examples of this process are shown in Fig.Q] Figure[2]shows 
all lines from all categories placed on the UVES intensity scale. 

(iv) Blend removal: If, when Gaussian-fitting a single line, an- 
other weaker blending line is present but ignored in the fit, then the 
centroid returned from the fit will be shifted towards the blending 
line. The magnitude of the shift will depend not only on the veloc- 
ity separation between the two lines, Av scp , but also on the relative 
intensities of the two lines, /2//1. When the two lines are not re- 
solved from each other, it is easy to estimate the velocity shift: the 
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Figure 2. Lines from all categories placed on the UVES intensity scale. 
The grey/red points are lines from the initial LPWN list which satisfy the 
blending criterion defined in step (iv) of the ThAr selection algorithm. The 
black points are lines satisfying all selection criteria and constitute the final 
ThAr line-list. Here, /]j st refers to the re-scaled listed intensities and 7 mca s is 
the amplitude of the line measured from the UVES ThAr spectrum. 



new centroid wavelength is approximately the intensity-weighted 
mean wavelength of the two blended lines, 



A c * 



+ l%Ai 



h +h 

The velocity shift due to the blending line is therefore 



A c - Ai 



sa Avv 



'I+/2//1 



(1) 



(2) 



However, if the two lines are further apart and are partially re- 
solved it is not clear how Av c will depend on Av sep and Fig- 
ure[3]rigfit) shows the results of a numerical experiment where two 
blended lines are varied in relative intensity and separation and fit- 
ted with a single Gaussian with similar initial guesses as in step 
(ii). A by-eye fit to the contours of constant velocity shift gives the 
following relationship: 



h 
h 



Av lol 



1 +0.1 



Av sf 



+ 0.0012 



Av si 



(3) 



The first term on the right-hand-side is just the previous equation 
and applies when Av sep is small with respect to <r v = (c/R)/2.355 
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Figure 3. Removal of strongly blended lines. Right: The shift in the centroid of a synthetic Gaussian due to blending with a weaker line. The black dashed lines 
are contours equally spaced in the logarithm of the shift. The solid grey/red line is a simple by-eye fit to these contours at a velocity tolerance of 40 m s" 1 and R = 
43000 [equation (3)]- All lines with blending lines more than 13 km/s away are safe from velocity shifts greater than 40 m s and this is marked with the vertical 
dashed grey/red line. This figure is available in colour in the online version of the journal on Synergy and from http://www.ast.cam.ac.uk/~mim/pub.html Left: 
Pairs of lines in the input line-list LPWN with the same solid and dashed grey/red lines from the right-hand panel. 



(c is the speed of light and c/R is the spectrograph's FWHM res- 
olution). We therefore reject lines from the input line-list LPWN 
which have blending lines within 13 km s -1 with relative intensities 
greater than that predicted by equation ([3} to produce shifts greater 
than a tolerance of 40ms~' at R = 43000. This tolerance is sim- 
ilar to the overall wavelength calibration residuals achievable at a 
typical UVES resolution of R = 43000 with the final ThAr list. 

(v) Removal of weak lines: If ThAr lines appear weak in the 
UVES lamp spectrum then the UVES pipeline's line-identification 
algorithm can fail or yield a false identification. The velocity pre- 
cision available in very weak lines is also too small to be useful 
in calibrating the UVES wavelength scale. We therefore reject any 
lines with measured intensities (above the measured continuum) 
less than 4 times the measured continuum level. 

(vi) FWHM selection: Any additional unknown features in the 
UVES ThAr spectrum can cause blending with the lines remaining 
after steps (iv) & (v) above. The next three steps aim to reject those 
lines which are effected in this way. The first of these steps is to 
remove lines whose widths are clearly inconsistent with the instru- 
mental resolution, in this case R as 70000 or FWHM as 4.3 km s~ . 
Figure [4] shows the distribution of FWHM for all lines which sur- 
vive steps (iv) & (v) above. After visual inspection of the lines lying 
away from the main cluster around FWHM~ 4.3 km s~ it was clear 
that lines wider than ~ 5.3kms~' or narrower than ~3.5kms~' 
should be removed. Lines were too wide when they were blended 
with other unknown features or where saturation of the CCD oc- 
curred. Lines were typically fitted as very narrow when they were 
very weak. 

(vii) High slope rejection: One of the improvements we made 
to the UVES pipeline was the addition of a continuum slope pa- 
rameter to the Gaussian fitting of ThAr lines during wavelength 
calibration. If a line is close to a strong, previously unknown fea- 
ture in the UVES ThAr spectrum then a large slope (relative to the 
line's intensity) might be needed to fit the line properly. Since these 
unknown features probably vary with lamp conditions, lines fitted 
with large continuum slopes are best avoided when calibrating the 
wavelength scale. Figure [5] shows the distribution of the absolute 
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Figure 4. The FWHM distribution of the lines surviving steps (iv) & (v) 
in the selection algorithm. Most lines lie in a well-defined band around the 
expected FWHM of ~4.3 ms"'. The dashed lines indicate our selection cri- 
teria for suitable calibration lines. The structure seen in this band is due to 
slightly differing resolutions in the different exposures from different wave- 
length settings. 



velocity difference, Av s , between two Gaussian fits - one made in- 
cluding the slope parameter, the other made with the slope fixed to 
zero - with continuum slope normalized by the line intensity. Vi- 
sual inspection of lines with Av s greater than 80ms~' showed that 
some of them have large asymmetries which are due to close blends 
with unknown features. Therefore, all lines with Av s > 80 m s~' are 
removed from the list at this stage. 

(viii) Large residual rejection: As a final selection step we re- 
ject lines which are fitted at wavelengths at some variance to those 
expected from the input line list LPWN. Since the ThAr line-list 
supplied with the UVES pipeline is known to contain several inac- 
curacies (see Section [TJ we do not apply this criterion in the first 
pass through the selection algorithm. In the second pass we reject 
lines which are fitted at positions more than 0.25 km s away from 
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Figure 5. Distribution of the absolute velocity difference between the cen- 
troids of two Gaussian fits - one made including the slope parameter, the 
other made with the slope fixed to zero - versus the normalized continuum 
slope. A conservative cut is made at Av 8 = 80ms~' to reject lines which 
might be effected by additional unknown features in the UVES ThAr spec- 
trum. 



Table 1. Excerpt from the final ThAr line-list for UVES. The complete 
table is available at http://www.ast.cam.ac.uk/~mim/pub.html and in the 
online version of this paper. The first column is the wavenumber (A vilc = 
1 X 10 8 /wo), the second column is t he air-wavele ngth computed from the 
vacuum wavelength, /l vac , using fhe|Edlen ( 1966) formula for th e disper- 
sion of air at 15°C and atmospheric pressure ( see iMurphv et al.ll200lL for 
detailed discussion). The third column is the logarithm of the original listed 
intensity normalized to the (arbitrary) intensity scale of the UVES ThAr 
spectrum. The fourth and fifth columns provide the line identification when 
available. The final column identifies the source of the wavenumber and 
which intensity scale the li ne wa s originally on: 'L' indica tes LP07, 'P' in- 
dicates[PE83, 'N' indicates |N73| while 'W indicates lw95l or lW02l Uniden- 
tified lines are generally given element designations 'XX' and ionization 
levels of '0' . For unidentified lines in 
if the line appears in PE83 or '1' if 
unknown. 



LP07 an ionization level '0' is given 



LP07 claim the line was previously 



co 


■tail 


login' 


Line Identity 


Source 


[cm" 1 ] 


[A] 




Element 


Ion 




26473.6290 


3776.271162 


3.299 


Th 


i 


P 


26450.5761 


3779.562436 


1.634 


XX 





P 


26426.4527 


3783.012697 


2.576 


Th 


ii 


P 


26424.4709 


3783.296425 


2.935 


Th 


ii 


P 


26415.5397 


3784.575599 


2.242 


XX 





L 


26408.3869 


3785.600691 


2.644 


Th 


II 


L 


26402.9308 


3786.383000 


2.653 


Ar 


II 


L 


26392.1220 


3787.933732 


1.813 


Th 


II 


L 



the expected position. In the third pass we reduce this parameter to 



0. 15 km s , less than l/10th of a pixel. 



4 FINAL THAR LINE-LIST 



The ThAr line-list formed after the third pass through the se- 
lection algorithm, referred to hereafter as the final or new list, 
is shown in part in Table [T] and is available in its entirety from 
http://www.ast.cam.ac.uk/~mim/pub.html and in the electronic ver- 
sion of this paper. In addition to an ASCII version, we also provide 



a version in a format which is readable by the UVES pipeline dis- 
tributed by ESO. 

The initial input list LPWN contains 13903 lines while the fi- 
nal list contains just 3070. For comparison, the line-list supplied 
with the UVES pipeline contains 2387 lines over the wavelength 
range covered by UVES. So, although a large number of lines are 
rejected in the selection algorithm, there are certainly enough re- 
maining lines to provide a reliable calibration of the UVES wave- 
length scale. Figure [6{left) shows the distribution of lines with 
wavelength in bins which approximate the size of extracted UVES 
echelle orders. Compared to the initial input list, the distribution 
of lines is quite uniform. This is mainly because of the rejection 
of close blends and demonstrates why the LPWN catalogue should 
not be used without first selecting unblended and otherwise reli- 
able lines. Of course, since LPWN contains so many lines, one can 
always achieve very low formal wavelength calibration residuals 
simply by defining a sufficiently small rejection threshold during 
the calibration. However, it is important to realise that with this 
approach many of the lines eventually used for calibration will ac- 
tually be blends which will distort the wavelength solution. 

Figure [6jleft) also shows that there are always more than 
~20 useful lines per UVES echelle order in the final list. This is 
enough to supply a robust wavelength calibration solution. Figure 
Upright) shows the distribution of lines at each step in the third 
pass through the selection algorithm. After the removal of close 
blends none of the subsequent steps removes lines in a strongly 
wavelength-dependent manner, as expected. Figure [TJleft) shows 
the contribution from Th, Ar and unidentified lines in the initial 
input list LPWN while Fig. |7jright) shows the composition of the 
final list. Note the strong decrease in the fraction of unidentified 
lines through the selection algorithm. This is bec ause m ost of the 
unidentified lines are those newly discovered by ILP071 many of 
which are quite weak and/or blended at the UVES resolution. 

We have used the final list to calibrate the UVES ThAr 
spectrum a final time and the resulting spectrum is available at 
http://www.ast.cam.ac.uk/~mim/pub.html and in the electronic ver- 
sion of this paper. In extracting and wavelength calibrating the 
spectrum with the modified UVES pipeline, we noted a dramatic 
improvement in the wavelength calibration residuals. To quantify 
this we wavelength calibrated the ThAr spectrum with two differ- 
ent inpu t line-lists - the final list derived above and the selection of 
DCH98 for R = 50000 — with the same pipeline parameters. The 
one exception was the tolerance level for rejecting discrepant lines 
in the pipeline's wavelength calibration step. This was adjusted to 
optimize the calibration for each list; i.e. a balance was sought be- 
tween the number of calibration lines and RMS residual per echelle 
order. The results are summarized in Fig. [8] which compares these 
two quantities for each echelle order in the 346- and 580-nm stan- 
dard UVES wavelength settings. 

Fig. [8] demonstrates that our new list provides a factor of >3 
improvement in the wavelength calibration residuals and, particu- 
larly for the 346-nm setting, provides more lines per echelle order 
on average. Note also the significant reduction in the number of 
orders - all at the bluest extreme of the 346-nm setting - which 
have less than 10 lines for calibration when using the new lis t. The 
346- and 580-nm settings also illustrate the impact of the new ILP07I 
ThAr catalogue on the residuals: the 346-nm setting co vers w ave- 
lengths not covered by HARPS (spectra from which the[LP07| cat- 
alogue was der ived) w hile the 580-nm setting uses lines almost ex- 
clusively from lLP07l The mean residual for the 346-nm setting is 
36 m s while that for the 580 setting is 28 ms^'-a small but clear 
improvement due to the improved line wavelengths in iLPOl Note 
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Figure 6. Histograms showing the expected number of lines per echelle order for different line-lists. Left: Comparison of the initial line-list with the list after 
close blends have been removed. Also shown is a comparison of the final line-list with that provided with the UVES pipeline. Right: The reduction of the 
line-list through the various selection stages. 
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Figure 7. Histograms showing the contributions from Th, Ar and unidentified lines to the initial (left panel) and final (right panel) ThAr line-lists. 
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Figure 8. Improvement of the wavelength calibration solution using the new 
line list (closed circles) over that achievable with the R = 50000 list from 
DCH98 (open triangles). The blue/darker points represent echelle orders 
from the 346-nm standard UVES wavelength setting while the red/lighter 
points are from the 580-nm setting. The calibration residuals for the 580-nm 
setting using our new list are slightly better than for the 346-nm setting due 
to the improved precision in the new |LP07| ThAr atlas. 



accuracies have on reduced spectra we have reduced the same raw 
ThAr exposures used above with the unmodified UVES pipeline 
and the ESO line-list. We then follow the line selection algorithm 
in Section [3T21 using a tolerance on the residuals of 0.3kms _1 in 
step (viii). Using the LPWN line-list as input allows us to trace 
any distortions of the wavelength scale introduced by the inaccura- 
cies in the ESO line-list and wavelength calibration procedure. The 
residuals between the Gaussian fit centroid wavelengths and the 
wavelengths listed in our final ThAr list derived above are plotted 
in Fig. [9] The mean residual shows strong, statistically significant 
variations with wavelength. The two most striking features are (i) 
the spike in the mean residual around 6300 A which has a peak-to- 
peak amplitude of ~75 m s and (ii) the large-scale decrease in the 
mean residual from ~0ms~' below 4500 A to — 30ms~' over the 
range 5500-7000 A and then the slow increase to — 1-30 m s -1 red- 
wards of 7000 A. These distortions of the wavelength scale have 
important consequences for several applications which have used 
pipeline-reduced UVES spectra. We illustrate one such application 
below. 



that the mean residuals in the two settin gs are the same (130m s -1 ) 
when the spectra are calibrated with the DCH98 list. 

Part of the motivation for this work was the various inaccura- 
cies we noticed in the ThAr line-list distributed by ESO with the 
UVES pipeline (see Section QJ. To investigate the effect these in- 



5 IMPLICATIONS FOR CURRENT UVES 

CONSTRAINTS ON A VARYING FINE-STRUCTURE 
CONSTANT 

Distortions of the wavelength scale of the magnitude and over the 
wavelength ranges just discussed are critically important to current 
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Figure 9. Residuals between the Gaussian fit centroid wavelengths mea- 
sured from the ThAr spectrum (reduced using the unmodified UVES 
pipeline and ESO line-list), Aa t , and the wavelengths listed in the LPWN 
line-list compiled in this work, /l|i st . Note the structure in the residuals as 
a function of wavelength. The mean residual, in bins which approximate 
the size of two UVES echelle orders, is shown as a grey/green thick line. 
The typical error in the mean residual in each bin is ~5— 10 ms" 1 . Thus, 
the strong variations from bin to bin and over longer wavelength ranges 
are statistically significant and indicate distortions of the wavelength scale 
introduced by the various inaccuracies in the ESO line-list and the UVES 
pipeline wavelength calibration procedures. 



UVES constraints on cosmological variations in the fine-structure 
constant, or em , which employ the UVES pipeline and ESO line-list. 
The velocity shift, v, of a transition due to a small relative varia- 
tion in a cm , Aa ja <sc 1, is determined by the ^-coefficient for that 
transition: 

r i n V Aa q 

co z = co + q \(a z /a ) 2 - 1 => - « -2 *- , (4) 

L J c a coq 

where ojq and co z are the rest-frequencies in the laboratory and at 
redshift z, a is the laboratory value of ff em and a- is the shifted 
value measured from an absorber at redshift z. The MM method 
is the comparison of measured velocity shifts from several tran- 
sitions (with different ^-coefficients) to compute the best-fitting 
Aa/a. Thus, the variation with wavelength in the magnitude and 
sign of the velocity shift in Fig.[9]will lead directly to spurious mea- 
surements of Aa/a depending on (i) which transitions are utilised 
to compute Aa/a and (ii) the redshift of the absorptio n system. 

An example of this problem arises in the results of Chan d^etall 

d2004h who studied 23 absorption systems and concentrated mainly 
on the Mgn /M2796/2803 doublet and the five strong Fen transi- 
tions redwards of 2340 A. The single absorber which dominates 
their sample, with a quoted uncertainty on Aa/a of just 0.1 x 10 5 , 
lies at redshift z abs = 1.2433. This places the three Fen transitions 
used in their analysis at observed wavelengths between 5250 & 
5340 A and the Mg n doublet appears at around 6270 A. Fig.[9]im- 
plies that a spurious velocity shift of 75 m s -1 between the Fe n and 
Mgn lines will be measured in this system. The three Fen lines 
all have large, positive ^-coefficients ~1400ctrr' while the Mgn 
doublet transitions both have very small (/-coefficients ~ 150 cm -1 . 
Equation Q therefore implies that the 75ms~' velocity shift cor- 
responds to a spurious (Aa/oOxhAr ~ +0.4 x 10~ 5 , four times larger 
than the 1-cr uncertainty quoted bv lChand et al.l d2004l) . 

This example of converting the residuals in Fig. [9] into pre- 
dictions of the systematic error in Aa/a are easily generalised to 
any absorption redshift once a set of n > 1 transitions is speci- 



fied. Fig. [10] shows two examples, one utilising the Mg n doublet 
and the five redder Fe n transitions, the other using those five Fe n 
lines plus the Feii/11608 transition. At each step in redshift, the 
value of (Aar/a^ThAr is calculated as follows: for each transition i, 
the 7 ThAr lines which fall closest to its redshifted wavelength are 
identified and their mean residual, v,, (as defined in Fig. [9j is cal- 
culated. The error in each v, is assumed to be the same and is set 
equal to the RMS residual of all 7x« ThAr lines, divided by V7. 
A linear least-squares fit between the v, values and the correspond- 
ing ^-coefficients [equation © then yields the best-fitted value of 
(Aa/Qr)xhAr and its uncertainty. 

The lower panel of Fig. QTJ sho ws the values of (Att /o-) ThAr 
for each of the 23 absorbers in the IChand et al.l ( 120041) analy- 
sis. As described above, some absorbers are particularly sensi- 
tive to miscalibrations due to their redshift and the particular tran- 
sitions used to determine Aa/a. However, for 18 of the 23 ab- 
sorbers, |(Aa/ar)ThArl is less than 0.2 x 10~ 5 . It is nevertheless im- 
portant to asses the impact this has on the overall weighted mean 
value of Aa/a. The values of Aa/a, corrected by their respec- 
tive (Att/a)xhAr values, are shown in the upper panel of Fig. [TTJ 
The errors bars are the quad rature addition of the original errors 
quoted bv lChand et alj ( 120041) and the errors in (Aor/o') ThAl shown 
in the lower panel. The weighted mean of the corrected values is 
Aa/a = (-0. 17±0.06)xl0~ 5 ; a2.5-cr 'detection' of variation in a em 
co mpared to the null re sult of Aa/a = (-0.06±0.06)x 10~ 5 claimed 
by IChand et al.l J2004I). Note that the error bars on the Aa/a val- 
ues from lChand et alj (2004) plotted i n Fig. 11 have been shown 
to be s ignificant underestimates in Murphy, Webb & Fla mbauml 
( 2007hf3). Thus, once larger, more appropriate error bars are used, 
the rel ative importance of wavelength calibr ation errors on the cor - 
rected lChand et alj fc004l) result diminishes dMurphv et al.l2007bl) . 
Nevertheless, Figs.[HjJ&[TT]clearly demonstrate how the ESO line- 
list and the wavelength calibration procedure of the original UVES 
pipeline cause significant systematic effects in measurements of 
Aa/a at the 0. 1 x 10~ 5 precision level. 

A similar analysis is presented by IChand et alj d2006h who 
cross-correlate two UVES ThAr spectra, one calibrated using the 
ESO line-list and the other calibrated using the line-list distributed 
by the National Optical Astronomy Observatory^ (NOAO), to de- 
termine the distortions of the wavelength scale over the range 
3100-6700 A. They find very similar structures (in their figure 7) 
to those seen in F ig. [9] including the sharp spike around 6300 A. 
IChandetal.ld2006l) then test the effect these distortions would have 
on a single Mg/Fen absorber at z-abs = 1.15, using only the red- 
der Fen lines and the Mgn double^. Since these lines are all 
redshifted to wavelengths where the mean calibration residual is 
similar (i.e. between 5040 & 6030 A) - falling just short of the 
spike around 6300 A - it is not surprising to find little-to-no ef- 
fect on Aa/a at th is redshift. Based on this single-redshift example, 
IChand et alj d2006t) claim that current UVES constraints on Aa/a 
should be fairly immune to any problems introduced by the ESO 



6 Although equation (4) formally has zero v-intercept, the fit between the 
Vj and q values must contain the y-intercept as a free parameter, to be de- 
termined simultaneously with (Aa/aOxhAr ■ This amounts to determining the 
best-fitted redshift of the set of the ThAr lines, mimicking the situation in 
the quasar spectra where one must determine the best-fitting Aa/a and red- 
shift simultaneously. 

7 http://www.noao.ed u/kpno/specatlas/thar/thar.html 

8 IChanaetalJfeOOd) do not consider the Fe n /41608 line in this test despite 
the fact that their analysis of the real absorber relies most heavily on this 
single transition. 
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Redshift 

Figure 10. The spurious value of Aa/a implied by the residuals in Fig. [9] (Aa/a)j^ r , as a function of redshift, assuming two different sets of transitions. The 
left panel utilises the Mg 11 ii2796/2803 doublet and Fe n i2344, 2374, 2382, 2586 & 2600 lines. The right panel uses the same Fe n lines plus the Fe n A 1608 
transition. The solid grey/green line shows a running median of (Aa/a;)niAr over 7 redshift bins. The redshift bins are Az = 0.01 wide. If one had measured 
a value of Aa/a from an absorber at a given redshift in a spectrum reduced with the unmodified UVES pipeline and calibrated using the ESO line-list, one 
would have to subtract the value of (Aa/a)jh Al at that redshift to derive a more reliable estimate. 
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Figure 11. Effect of wavelength calibration errors on the Aa/a values of 
IChand et alj J20041) . The lower panel gives the systematic error (Aff/aOThAr 
for each absorber while the upper panel shows the original Aa/a val- 
ues (red/grey squares) corrected by their respective (Aa/a^Ar values 
(solid circles). Even though the weighted mean calibration error is just 
(Aff/oOThAr = (+0.02 ± 0.02) X 10~ 5 , the fact that some of the (apparently) 
most precise Aa/a values have large corrections means that the corrected 
weighted mean Aa/a shifts significantly to Aa/a = (-0.17 ± 0.06) X 10~ 5 . 



line-list and (unmodified) UVES pipeline wavelength calibration 
procedures. As demonstrated in the specific examples above and in 
the more general cases illustrated in Figs.[T0]&[TT] such a simplis- 
tic extrapolation to other redshifts significantly underestimates the 
systematic errors involved. 

Finally, note that in iMurphv et alj fcOOll) and [Murphy et al. 
J2003I) we carried out a similar analysis on the Keck/HIRES spectra 
which provided our previous constraints on Aa/a. The wavelength 
calibration employed the NOAO ThAr line-list discussed above. In 
those works we found no such structures like those seen in Fig. 1101 
Thus, although the current UVES constraints on Aa/a have been 
adversely affected by the ESO line-list, such a problem does not 
generally exist for HIRES constraints to any significant degree. 



6 CONCLUSIONS 

We have implemented a simple iterative algorithm to select ThAr 
lines for use in calibrating echelle spectrographs. The algorithm 
combines the expected line properties (e.g. degree of blending), 
based on the input line-list, with properties of the ThAr spectrum 
observed on the spectrograph of interest (e.g. resolution, relative 
line intensities etc.). We emphasise that although this paper has fo- 
cused on the application of this algorithm to VLT/UVES spectra, 
the algorithm can be applied to any other spectrograph; all that is 
required is a ThAr spectrum taken on that spectrograph and a re- 
duction code which allows repeatable wavelength calibration with 
different ThAr line-lists. 

Using the most rec ent and reliab le ThAr laboratory wave- 
lengths, including that o f lLovis & Pepel J2007t) . we have produced 
a new ThAr line-list which, we argue, should be used to calibrate 
UVES spectra in preference to the list distributed by ESO. The 
main advantages of the former over the latter are: (i) more homo- 
geneous input line catalogues with mutually consistent wavelength 
scales; (ii) more accurate and precise laboratory wavelengths; (iii) 
only unblended lines appear in the final list; and (iv) only lines 
which are observed to be useful in UVES spectra are included. In- 
deed, the wavelength calibration residuals achieved using the new 
UVES line-list (RMS ~ 35 m S" 1 ) are mo re than a factor of 3 bet- 
ter tha n those achieved with the ESO or lDe Cuyper & Hens berge 
(1998) list (RMS ~ 130 ms" 1 ) - see Fig. [8] An assumption one 
makes in using our new line-list to calibrate UVES spectra is that 
the relative line intensities (i.e. lamp conditions like pressure, cur- 
rent etc.) are similar to those in the sample UVES ThAr spectra 
used in our selection algorithm. However, ideally, the selection al- 
gorithm should be included in the wavelength calibration procedure 
itself, thus producing the best selection of ThAr lines for each indi- 
vidual ThAr exposure being reduced. 

The new line-list has been used to probe distortions of the 
wavelength scale introduced into UVES spectra which were orig- 
inally wavelength calibrated with the ESO line-list. The distor- 
tions include sharp, short-range corrugations with peak-to-peak 
amplitude up to ~75ms~' and longer-range variations between 
+30ms~' and -30ms~' (Fig. [9j. These distortions have serious 
implications in the context of current UVES estimates of possi- 
ble variations in t he fine-structure constant, especially those of 
IChand et al.1 d2004l) - see Figs.[l0]&[TT] For example, if one em- 
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ploys the common Mg/Fen combination of transitions to derive 
Aa/a then the errors in the ESO line-list cause a spurious value 
of (Aa/aOxhAr ~ +0.4 x 10~ 5 for absorption redshifts 1.22 <z a b s ;S 
1.35 a nd (Att/«) ThA r ~ - 0-3 x 10~ 5 for 1.38 < z abs < 1.75. Al- 
though Chand et al. (2004) claim an overall precisi on in Aa/a of 
0.06 x 10~ 5 from 23 Mg/Fen absorbers (but see iMurphv et al.1 
l2007bl l3). a single absorber at z abs = 1.2433 strongly dominates 
their sample with an individual error on Aa/a of 0.1 x 10~ 5 . How- 
ever, at this redshift, and considering which transitions are used to 
derive Aa/a in this system, we find that a spurious systematic shift 
of (Aff/oOThAr ~ +0.4 x 10~ 5 has affected this system. Once such 
wavelength calibration corrections are applied to the entire sam- 
ple, the null result of Aa/a = (-0.06 ± 0.06) x 10~ 5 claimed by 
IChand et al.1 j2004l) becomes Aa/a = (-0. 17 ± 0.06) x 10~ 5 - a 2.5- 
cr 'detection'. The new ThAr line-list we present here should aid in 
avoiding such important yet elementary errors in future analyses of 
UVES spectra. 
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The following supplementary material is avail- 
able for this article online on Synergy and at 
http://www.ast.cam.ac.uk/~mim/pub.html 

Table [U The final ThAr line-list for UVES. In addition to an 
ASCII version, we also provide a version in a format which is read- 
able by the UVES pipeline distributed by ESO. 
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